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The phase diagram of the unconstrained t — J model is calculated using the 
random phase approximation. It is found that the extended s and the d x 2_ y 2 -channels 
are not degenerate near half filling. Extended s-pairing with a low T c occurs only 
for a band containing less then 0.4 electrons or holes per unit cell, whereas in a 
large region around half-filling d-wave pairing is the only stable superconducting 
solution. At half filling superconductivity is suppressed due to the formation of the 
anti-ferromagnetic Mott-Hubbard insulating state. By extending the analysis to the 
unconstrained t — t' — J model, it is proven that, if a Fermi surface is assumed 
similar to the one that is known to exist in cuprous oxide superconductors, the 
highest superconducting T c is reached for about 0.7 electron per site, whereas the 
anti-ferromagnetic solution still occurs for 1 electron per site. It is shown, that 
the maximum d-wave superconducting mean field transition temperature is half the 
maximum value that the Neel temperature can have in the Mott-insulating state. 
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In spite of a huge experimental and theoretical effort to understand the superconductivity 
in the cuprous oxide high T c superconductors []1],|2|, a mechanism for superconductivity has 
not yet been firmly established. A lot of attention has been devoted to obtain a model 
for superconductivity starting from the Hubbbard hamiltonian, however, there is a growing 
suspicion that the positive U Hubbard model alone can not give rise to superconductivity 
||. A different approach has been to treat the electronic degrees of freedom and 'spin 
fluctuations effectively as separate channels leading to a retarded electron-electron 

interaction mediated by spin-fluctuations. The latter model has proven to be more successful 
in the sense of providing a possible mechanism of superconductivity. 
Here I will discuss pairing using an exchange-only Hamiltonian 

+ 2 J^Q J2k,qJ(Q) Sk,k+Q ■ Sq,q-Q ~ \ J2aa' C k+Qa C la' C q+Q^ c ka 

The t — J model is studied here without the usual constraint on double occupancy of 
the same site as a model in its own right. In the real-space representation J couples 
nearest neighbouring sites on a square lattice. Hence the exchange part is of the form 
J2r,s <3r ■ Sr + $ — \n R riR + s) where 5 runs over the four nearest neighbour sites, and the 
factor 1/2 compensates double counting of the interactions in the summation over R. If 
one considers only two neighbouring spins the energy of the triplet state is 0, and of the 
singlet it is — Jo, hence the sign convention is such, that J > couples neighbouring spins 
anti-ferromagnetically. 

For the dispersion of the band I will consider a nearest neighbour and a next nearest neigh- 
bour hopping term, but apart from that, it is the same hamiltonian as was used in the 
papers by Baskaran, Zou and Anderson (BZ A) |||7| , Emery and by Kotliar . Although 
usually such a Hamiltonian is thought to be derived from the Hubbard model by means 
of a Gutzwiller projection, (which also changes the character of the fermion operators, by 
projecting out double occupancy of the same sites) it should be pointed out, that in the 
cuprous oxide systems this term may also have a different microscopic origin. As the ac- 
tual bandstructure in these systems is experimentally known to be better described by the 



three band model of Zaanen, Sawatzky and Allen []10| , (which is again a simplified version 
of the real valence band structure involving 6 oxygen 2p bands and 5 cupper 3d-bands for 
the occupied states, as well as unoccupied 3s and 3p states) a transformation to a single 
band hamiltonian will in principle generate both an effective Hubbard U and an intersite 
J pT^T5[. Examples of such transformations can be found in the work by Emery |J, and 
by Jansen jT4|. However, also other, more complicated types of interactions are generated 
when making transformations of this type, notably the correlated hopping term (with six 
operators) which, as has been shown by Hirsch, promotes superconductivity of hole-carriers 
|T5| . The interaction considered by Jansen as well as the correlated hopping term treated 
by Hirsch, effectively provide an on-site attraction, which, when considered on its own, 
promotes pairing in the (non-extended) s-wave channel. In this paper I will assume that 
the net on-site interaction is repulsive, which, as will be discussed below, tends to suppress 
superconductivity by stabilizing the anti-ferromagnetic solution. 

Monthoux and Pines have considered the t—t' bandstructure together with an interaction 
of the form g{q)s(q) ■ S(—q), where s(q) represents the valence electron spin-operator, and 
S(—q) is a separate spin-fluctuation operator, the properties of which are determined by the 
spin-susceptibility. The transport and superconducting properties are then calculated from 
strong-coupling theory using empirical values for the spin-susceptibility and g(q). It has 
been shown by Monthoux el al. that the effective electron-electron Kernel arising from such 



a coupling becomes indeed a spin-dependent interaction (T^|, which could in principle be 
written as a frequency dependent version of Eq. p] [|17j] . In the same paper a weak coupling 
analysis of such a frequency dependent Kernel was given. In the present paper the prob- 
lem is further simplified by disregarding a possible frequency dependency of J(Q), which, 
depending on the details of the microscopic origine of Eq. |1], may still be a justifyable ap- 
proximation. This allows us to explore the phase diagram in somewhat more detail without 
having too many parameters to consider. 

Apart from these general considerations I have no sound microscopic justification for using 
this hamiltonian. The main motivation to use it comes from the fact that, as I will show in 



this paper, it appears to do a surprisingly good job as a phenomenological model consistent 
with at least some of the experimentally known aspects of superconductivity in these sys- 
tems. 

BZA |7j considered pairing of the s*-type near half filling, Emery considered d x % -^-pairing, 
and Kotliar studied both s* and d-type pairing. Below I will show, that the s*-type pairing 
is not a stable solution near half filling, and is dominated by pairing of the cf-type. As the 
latter again tends to be unstable with respect to the anti-ferromagnetic Mott-Hubbard in- 
sulating state at half filling, superconductivity can only exist sufficiently far away from this 
region. As the optimal T c would have been reached at half filling for a symmetrical band, 
this would lead to the conclusion that superconductivity is only a marginal effect in such 
a system. However, the high T c cuprates do not have a symmetrical band, and the Fermi 
surface is known to be distorted from the perfect square that arises from considering only 
nearest neighbour hopping. This actually comes to rescue: As a function of band- filling it 
pulls apart the regions, where anti-ferromagnetism and high T c are optimal, without hav- 
ing a noticable effect on the superconducting or anti-ferromagnetic transition temperatures. 
Three important trends emerge from this analysis: 

(1) Given the distorted shape of the Fermi surface as it is known to occur in the cuprates, 
hole doping gives rise to higher T c 's than electron doping. 

(2) There exists a universal relation between the highest Neel temperature found in the 
phase diagram and the highest possible mean field superconducting T c , with T N /T C « 2. 
That a relation of this kind should exist was already pointed out by Anderson [[SJ shortly 
after the discovery by Bednorz and Miiller. 

(3) This implies that with reasonable values for the intersite exchange interaction, provid- 
ing the correct Neel temperature, one automatically obtains values for the superconducting 
transition temperature which are (even though they are larger than the experimental values 
in the cuprates) definitely in the right ball-park. 

The fc-space representation of the exhange interaction is of the form of Eq.|l| with 
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J{Q) = 2 J o( cos (Qxd) + cos(Q y a)) (2) 

This type of interaction favours anti-ferromagnetism if Jo > 0, which becomes especially 
stable if the band is half filled. The antiferromagnetic alignment of nearest neighbours gives 
rise to a spin-dependent effective field, which is periodic with the wave vector (ir/a, 7r/a). 
Let us now turn to the superconducting gap equations. If the interaction potential Vk q is of 
the form 



H l = ~ V k,q°h C -ki C -<li C <l\ ( 3 ) 
k,q 



the BCS gap equation is |24| 



Y 2E i 2ksT v ; 

where E q = ^/e|~+~A| as usual. If we can make the assumption, that the main contribution 
leading to superconductivity comes from the J-term, we see that the interaction entering 
the gap equations is 

V k>q = 2J(k -q) + 2J(k + q) (5) 
With this substitution we obtain 

Afe = J ^ cos k x a + cos k y a (6) 
where I introduced the dimensionless pairing amplitudes 

= E g ^ 1 A 9 cosg i atanh( 1 g T ) (7) 

As there are two possible order parameters tyf, we have here two coupled equations, which 
can be easily disentangled with the help of symmetry selection rules. I will do this for the 
case where the superconductor has a four-fold rotation axis. In that case A is either an odd 
or an even function of k. In the former case, which corresponds to d xy symmetry, pairing 
amplitudes of the form = J2 q Eq 1 ^ sin q^a have a finite amplitude, whereas and 



are zero. The hamiltonian considered here does not couple to the d xy pairing-channel. If A 
is even, we have to consider two possibilities: Either = leading to an extended s-wave 
gap, or = — tyy leading to a d x 2_ y 2 symmetry of the gap function. The gap function 
corresponding to these two cases is 



Afc = -AofcosA^a ± cos k y a] (8) 



where the plus and minus sign correspond to the s*- and J^-^-wave types of pairing re- 
spectively, and the gap equation becomes 

2 E 

— = ^E- 1 [cosq x a±cosq y a] 2 t&nh(^-^) (9) 

This equation was also obtained by Kotliar ||. In his analysis the constraint of no 
double occupancy of the same site was taken into account in an appriximate way, by having 
t proportional to doping of the half filled band. At half filling one then effectively has t = 0, 
for which case, as was shown by Kotliar, the summations on the right hand side of this Eq. 
|9] are identical for the two types of pairing. As a result he obtained a degeneracy between 
the d-and s-ordered state at half filling, leading to the conclusion that a pairing of type 
s + id could occur. For any finite value of t this degeneracy is however lifted. In the mean 
time a variety of numerical and theoretical techniques have been applied to the t — J and 
related models, from which a tendency toward d x 2_ y 2 -pairing has been found near half filling 
25| f2~9|. It is easy to show, that for the s*-type pairing at half filling of a symmetric band, 



Jo has to exceed a critical value. Let us assume that e = (W/4) ■ (cos(k x a) + cos(k y a)). For 
T = the gap equation becomes 



J 1 ^/W 2 /A + Al = \cos q x a + cos q y a\ = 0.811 (10) 
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hence the cricial value of Jo for s*-pairing is at half filling Jq = 0.62W. The reason for the 
appearance of a critical value is, that at half filling the s*-type gap is exactly zero for all k at 
the Fermi surface. Only by mixing in states away from the Fermi level, superconductivity of 
this type may occur, which requires a minimum coupling strength. The <i-channel is much 
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more effective in this sense, as is finite at the Fermi surface except for the node-points. 



As a result in the (i-channel we have Jq = 0. 



The ground state energy relative to the normal state can now be determined by realizing 
that it is the expectation value of the reduced hamiltonian minus the non-interacting part, 
which is |St| 



where the first two terms represent the energy gained by redistributing the electrons over 
/c-space in the correlated wavefunction, whereas the third term compensates double counting 
of the interaction. In principle one has to solve the gap equation together with a constraint 



on electron occupation number ||3T|-|33|j, however the corrections to the free energy are of the 



order (A /Ef) 2 , [|34] which is small for the parameters that we will consider. 



I still need to specify the electron dispersion relation before we can solve the gap equations. 
If one considers a tight-binding model with a single orbital per site, with only hopping 
between nearest and next nearest neighbours, the single particle energies are 

efc = —It (cos(k x a) + cos(k y a)) — At' cos(k x a) ■ cos(k y a) — E F (12) 

The t'-term is due to next-nearest neighbour hopping. Let me briefly discuss some of the 
properties of such a band. If t' — at half filling of the band, such a dispersion relation 
has the remarkable property that the Fermi surface forms a perfect square, with a diverging 
effective mass over the entire Fermi surface. In practice this situation will never occur, as 
there will always be some finite coupling between next nearest neighbours. This causes a 
bulging of the Fermi surface, as is shown in Fig. [j], which eventually transforms into a 
rotated Fermi surface if \t'\ 3> \t\. The shape obtained for t' = —0.7t is very close to what 
has been calculated with the local density approximation for e.g. La2Cu04 and YBa 2 Cu307 



35| , |36[| . A significant change also occurs in the density of states (DOS) at the Fermi energy, 
which is displayed in Fig. ^ as function of the number of electrons per unit cell. This 
somewhat unusual representation of the DOS is useful in the discussion below, where we 



compare ground state energies of various types of ordering at a given electron density. We 
see, that as t' is increased, the DOS becomes a-symmetric, and the maximum is shifted to 
the left side of the point where the band is half filled. Of course the direction in which this 
occurs is dictated by the sign of With t' < we mimic the situation encountered in the 
CwCVplanes of the high T c cuprates. 

In Fig. § numerical calculations of the free energy are shown as a function of occupation 
number for Jq/W = 0.6, where W = 8t corresponds to the bandwidth if t' = 0. For 
the sake of completeness also the free energy of the anti-ferromagnetically ordered state is 
included. This was calculated from the same hamiltonian. To stay in the same spirit as 
for the superconducting solutions, the random phase approximation was used. Hence the 
free energy was minimized together with a constraint on the electron occupation number, 
anticipating a finite expectation value of < c\ + q^c^ >= — < c -ki c -k-Qi > a t the point 
Q = (7r/a,7r/a). We notice that the anti-ferromagnetically ordered state at half filling is 
always more stable than the metallic state. However, for small values of Jo the ci-wave 
paired superconducting state is still more favourable. This is a consequence of our choice of 
model Hamiltonian, which is perhaps somewhat pathological near half filling: Physically the 
exchange terms should arise from a strong repulsive interaction between electrons making 
a virtual transition to the same orbital of e.g. a transition metal atom. On the one hand 
this leads to exchange coupling between spins on neighbouring orbitals, while on the other 
hand it causes the opening of a Mott-Hubbard gap, which is much larger than the anti- 
ferromagnetic gap. This would strongly stabilize the anti-ferromagnetic solution. Tempting 
as it may be to add an on-site repulsion at this point as an additional model parameter, I 
will not do so: It has become clear in recent years, that a large repulsive U gives rise to 
very strong correlation effects, at and near half filling, which can not be properly treated 
with the random phase approximations made in this paper []37| p9f| ■ For this reason, and 
also because fluctuations are neglected with the latter approximations, the present analysis 
is insufficient close to half filling. For higher doping it could have some relevance to the 
mechanism of superconductivity. It is important to add in this context, that the symmetry 
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of an additional on-site interaction is such, that it cancels out in the gap equation for the s* 
and rf-channels. Hence an on-site U does not affect the gap-function or the free energy for 
these types of superconductivity 

Although from a Maxwell construction one is lead to the conclusion that phase separation 
should occur in s*- and d-ordered regions, this is strongly suppressed if the long range 
Coulomb interaction is taken into account. [18] Although the Coulomb term is not included 



explicitly in the Hamiltonian, the presence of such a term is assumed implicitly by imposing 
the constraint that the electronic density is macroscopically conserved. As was stressed by 



Emery, Kivelson and Lin flHy^]], who studied the t-J model together with the constraint 
on no doubly occupied sites, 'the holes are often donated by oxygen atoms which are quite 
mobile ... ', providing a physical mechanism for screening of the long range Coulomb term. 



Putikka, Luchini and Rice |21| provided numerical evidence that, in the absence of a long 



range Coulomb force, phase separation occurs for J/t > 3.8 as n — > and J/t > 1.2 near 
half filling. The present analysis does not lead to phase separation if the long range Coulomb 
interaction is taken into account. 



The phase diagram is displayed in Fig. Due to electron-hole symmetry in this case, the 
diagram is symmetric around half occupation of the band. Roughly speaking s*-pairing is 
favoured far away from half filling of the band, whereas c?-wave pairing becomes the most 
stable solution near half filling. For J < 0.31^ there are regions of no superconductivity, 
which broaden upon decreasing J , and completely cover the horizonal axis for precisely Jo — 
0. This tendency towards <i-wave pairing near half filling was also obtained by Littlewood 
40| for the charge transfer model | I0| , again using a weak coupling treatment. In these 



calculations an inter-site exchange interaction is not introduced explicitly, and can only 
result inderectly from the repulsive on- and inter-site interactions which are taken into 
account in the model. 

We see, that Jo > 0.7W is required to find an antiferromagnetic phase near half filling. As 
can be seen from the free energies versus doping, the phase boundaries between s* and d, and 
between d and AF, correspond to a discontinuous change from one type of ordering to the 
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other. For the s*-d boundary this discontinuity will probably be softened without loosing 
the superconductivity by the occurrance of an intermediate state of mixed s + id character, 
as was proposed by Kotliar at precise half filling and t — 0. The phase boundary between 
d and AF is different in this respect. As both the anti- ferromagnetic and superconducting 
correlations occur in the same band of electrons, they will tend to suppress each other. 
Because finite anti-ferromagnetic correlations will occur on the superconducting side of the 
phase boundary and vice versa, at the boundary T c and T/v should come out to be zero if 
such corrections are taken into account. This requires a treatment of the model hamiltonian 
which goes beyond the level of random phase approximations made in this paper. The fact, 
that the (i-paired and anti-ferromagnetic solutions both have their optimum at half filling, 
is rather worrying, as in a real solid the anti-ferromagnetic solution will in practice turn out 
to be the more stable one, due to the opening of a Mott-Hubbard gap. 
Fortunately nature does provide us with a way to make a separation in parameter space 
between the anti-ferromagnetic and superconducting states. As already pointed out above, 
in practice there will always be a finite value of t'. From Fig. Q we see, that in this case the 
maximum value of the DOS does not occur at half filling of the band. A well-known result 
from BCS theory is, that a high DOS at the Fermi level enhances T c . If on the other hand 
precisely 50 % of the states is occupied, the opening of an anti-ferromagnetic gap causes 
a downward shift of all occupied levels in the reduced Brillouin-zone, which is the reason 
why the anti-ferromagnetic solution is best stabilized at precise half filling. This effect is 
demonstrated in Fig. [5], where we see that indeed the lowest free energy of the <i-paired state 
occurs now at 35% filling, whereas the anti-ferromagnetic solution is still at half filling. We 
also notice from this plot, that this a-symmetry implies that the highest T c 's of a <i-paired 
superconductor are to be expected on the left ('hole-doped') side of half-filling. Lower T c 's 
occur on the right side. 

Let us now consider the A/T c -ratio following from the gap equation. Within the context of 
BCS theory we have A (T) = at T c , so that T c follows from 
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— = ^e g 1 tanh( q )[cos q x a ± cos q y a] 2 (13) 

^0 q *KbJ-c 

where the ± sign refers again to the two symmetries of pairing. This equation can be easily 
solved numerically The result is, that for extended s-wave pairing the ratio 2Ao//cbT c is 
6.5, whereas for d-wave pairing it rises gradually from 4 if J <C W, up to 6.5 in the limit 
where Jo 3> W. This is not sensitive to the value of the parameter t' . We should keep in 
mind here, that Ao is the maximum value reached by A (A;) (respectively at the (it, 0)- and 
(7r, 7r)-point for d- and s*-pairing). 

Finally it is interesting to look how the mean field estimate of T c depends on the coupling 
strength J /W. In Fig. |T C MF /W is displayed function of Jq/W for the d-wave channel. 
First of all we notice, that for J > W/4 the value of T C MF is about J /4. For J /W « 1 
this crosses over to a quadratic dependency T C MF = AJq/W . For comparison a similar curve 
is displayed for conventional s-wave pairing, using the negative U Hubbard model in a band 
with a square DOS. We notice that the mean field transition temperature with the latter 
model becomes T C MF = \U\/A for large \U\ (which is actually outside the range of validity 
of the BCS weak coupling approach |4Ti^2l ), and has the familiar BCS-like exp (—W/\U\) 
behaviour for small U . The T c for the extended s-wave pairing lies again below the negative 
U curve, and is only finite above a threshold value of Jo as discussed above. 
A consequence of this is, that the model hamiltonian proposed here leads to quite reasonable 
values of the transition temperature, which are relatively insensitive to the value assumed 
for the bandwidth. If we assume that for example J = O.leV, we would find that the Neel 
temperature in the Mott-insulating state can not exceed the mean field value for Z = 4 
interacting neighbours T^ IF = ZS(S + l)Jo/6 = 580K. If we assume that the bandwidth is 
smaller than about 0.5 eV, we obtain from the BCS gap equation that the d-wave transition 
temperature can not exceed the mean field value T C MF = 290K. This demonstrates that the 
optimal Neel temperature and the optimal superconducting transition temperature have a 
ratio of about 2 for Jq/W of the order 0.2 to 1. Keeping in mind, that with the mean field 
approach we over-estimate both and T c , I expect that the ratio between the two should 
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remain relatively intact if corrections beyond the mean field approximation are included. 
In the limit where Jq/W is small, T c comes out smaller, although the suppression of the 
transition temperature goes much slower then for conventional s-wave superconductivity. 
For example if the bandwidth is 1 eV we obtain T c = 137 K, and with W = 2 eV we find 
that T c = 74 K. 

Finally it is possible now to draw a phase diagram in the temperature versus doping 
plane. Let us choose Jo = O.leV, which gives approximately the correct value for Tn, 
and t' jt = —0.7 which gives approximately the right Fermi surface. To stay in the regime 
where T/v ~ 2T C , let us assume 8t = W = 0.5eV. The latter parameter is rather small 



compared to the 2 to 3 eV of the bare copper-oxygen p^, p y , d x 2_ y 2 anti-bonding band j35 
and leads to a slight over-estimation of T c . Although the estimated Neel temperature can 
be indicated at half filling of the phase-diagram, the present analysis has no bearing on the 
region near half filling, which was left open for that reason. The phase diagram, displayed 
in Fig.|7|, perhaps somewhat optimistically gives values of T c above 200 K, which sofar has 



not been found experimentally. According to the rule-of-thumb that T c scales with T/v P3 



one has to look for systems with relatively high Neel temperatures in order to reach room 
temperature superconductivity. A number of factors will lower T c below the mean field value 
given here. First of all, anti-ferromagnetic correlations will occur near the Mott-insulating 
state, which tend to suppress the superconducting order. Very strong on-site spin corre- 
lations are known to exist due to the large on-site U, but d and s*-pairs are insensitive 
to this interaction channel, as can be seen from the gap equation. Second, the mean field 
approach links T c directly to the energy scale of the pair-breaking, which again leads to an 
overestimation of T c . The reason for this, is that the long range phase- coherence can be 
lost in a dephasing-transition, if the phase fluctuations have a lower energy scale than the 
pair-breaking energy. This requires a better knowledge of the phase fluctuation spectrum of 
the d-wave superconducting state. 

Using a weak coupling BCS treatment of the t-t'-J model, I have shown that there exists a 
universal ratio of 2 between the Neel temperature at half filling and the optimal mean field 
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superconducting transition temperature. If a realistic shape of the Fermi surface is taken, 
the optimal T c occurs for 0.7 electron per site, while the Mott-insulating antiferromagnetic 
state occurs at half filling. With these parameters, T c is shown to be lower for electron 
doping than for hole doping. 

Acknowlegments It is a pleasure to thank Z. X. Shen for a stimulating discussion and com- 
munication of some unpublished results, and J. Lorenzano and G. A. Sawatzky for useful 
comments on the manuscript. This work is part of the research program of the "Stichting 
voor Fundamenteel Onderzoek der Materie", which is financially supported by the "Neder- 
landse Organisatie voor Wetenschappelijk Onderzoek". 



13 



REFERENCES 
[1] J. G. Bednorz and K. A. Miiller, Z. Phys. B. 64,189 (1986). 
[2] M. K. Wu et al , Phys. Rev. Lett. 58, 908 (1987). 

[3] A. Sudbo, C. M. Varma, and T. Giamarchi, Phys. Rev. Lett. 70, 978 (1992). 
[4] D. J. Scalapino, E. Loh, and J. E. Hirsch, Phys. Rev. B 35, 6694 (1987). 
[5] P. Monthoux, and D. Pines, Phys. Rev. B 47, 6069 (1993). 

[6] P. W. Anderson, G. Baskaran, Z. Zou, and T. Hsu, Phys.Rev. Lett. 58, 2790 (1987). 

[7] G. Baskaran, Z. Zou, and P. W. Anderson, Solid State Commun. 63, 973 (1987). 

[8] V. Emery, Phys. Rev. Lett. 58, 2794 (1987). 

[9] G. Kotliar, Phys. Rev. B 37, 3664 (1988). 
[10] J. Zaanen, G. A. Sawatzky, J. W. Allen, Phys. Rev. Lett. 55, 418 (1985). 
[11] F. C. Zhang, and T. M. Rice, Phys. Rev. B 37, 3759 (1988). 
[12] H. Eskes, and G. A. Sawatzky, Phys. Rev. Lett 61, 1415 (1988). 
[13] J. H. Jefferson, H. Eskes, and L. F. Feiner, Phys. Rev. B45, 7959 (1992). 
[14] L. Jansen, Physica C 156, 501 (1988). 
[15] J. E. Hirsch, Phys. Rev. B 43, 11400 (1991). 

[16] P. Monthoux, A. V. Balatsky, and D. Pines, Phys. Rev. B46, 14803 (1992). 

[17] A small difference is the absence of the spin- independent second term of Eq. [I] in the 
expressions of Monthoux. It implies that in their analysis the triplet channel is repulsive, 
and the interaction in the s* and d x 2_ y 2 channels is scaled with a factor 3/4 compared 
to our expressions. Also a repulsive interaction in the d^-channel emerges. 

[18] The charging energy per unit volume for a small fraction of spherical d-type droplets 

14 



dispersed in an s*-type medium is j^Tx 2 e 2 (n d — n) 2 r 5 • N/V, where N = " d ~" 3 4 J^/ 3 is 
the number of droplets, V is the sample volume, rid, and n s are the densities of the d 
and s-type ordered states, and n is the average density. For n d and n s of comparible 
magnitude we interprete r as the characteristic length scale on which phase separation 
occurs with a charging energy per unit volume 2 fe 2 {rid — n)(n — n s )r 2 , which has the 
correct limiting behaviour for small n s or n d and interpolates between these two cases. In 
a layered electron gas one can instead consider a cilindrical shape, for which the charging 
energy density is ^e 2 (n d — n)(n — n s )r 2 . The r 2 scaling behaviour tends to reduce r. 
A lower bound occurs due to quantum size effects: Due to quantum confinement the 
kinetic energy of each electron increases with an amount of about n \^lf ■ For particles 
moving in a 2D plane Ak 2 ~ so that the increase in energy per unit volume is 
(nd — n)(n — n s )(rid + n s — n)(n^ — n s )~ 2 , which again interpolates between the two 
limits where or n s is small. 

[19] V. J. Emery, S. A. Kivelson, and H. Q. Lin, Phys. Rev. Lett. 64, 475 (1990). 

[20] S. A. Kivelson, V. J. Emery, and H. Q. Lin, Phys. Rev. B 42, 6523 (1990). 

[21] W. O. Putikka, M. U. Luchini, and T. M. Rice, Phys. Rev. Lett. 68, 538 (1992). 

[22] The minimum energy due to the combined effect of the Coulomb interaction and quan- 
tum confinement occurs at a length scale r 4 = 8ao(nd + n s —n){nd — n s )~ 2 . The increase 
in energy per unit volume is then 7rRy(nrf — n){n — n s )(rid — n^ 1 ^ (2a ) 3 (nd + n s — n), 
where Ry and a are the effective Rydberg and Bohr radius respectively. This should 
be smaller than the gain in free energy due to phase separation. For the parameters 
considered in this text, the latter is at least two orders of magnitude smaller than the 
former, implying that phase separation is strongly suppressed. 

[23] P. W. Anderson, science 235, 1196 (1987). 

[24] J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Phys. Rev. 108, 1175 (1957). 



15 



[25] C. Gros, R. Joynt, and T. M. Rice, Z. Phys. B 68, 425 (1987). 

[26] G. Kotliar, and J. Liu, Phys. Rev. B 38, 5142 (1988). 

[27] M. Grilli, and G. Kotliar, Phys. Rev. Lett. 64 1170 (1990). 

[28] M. Grilli, R. Raimondi, C. Castellani, C. di Castro, and G. Kotliar, Phys. Rev. Lett. 
67, 259 (1991). 

[29] E. Dagotto, and J. Riera, Phys. Rev. Lett. 70, 682 (1993). 

[30] M. Tinkham, 'Introduction to superconductivity', Mc Graw Hill (New York) (1975). 
[31] G. Rietveld, N. Y. Chen, D. van der Marel, Phys. Rev. Lett. 69, 2579 (1992). 
[32] D. van der Marel, Physica C 135, 35 (1990). 

[33] D. van der Marel, and G. Rietveld; Phys. Rev. Lett. 69, 2575 (1992). 

[34] D. van der Marel and G. Rietveld, in 'Electronic properties of high T c materials. The 
normal and superconducting state of High Tc materials.', Editors: H. Kuzmany, M. 
Mehring, J. Fink, Springer Series in Solid State Sciences, (1993). 

[35] W. E. Pickett, Rev. Mod. Phys. 61, 433 (1989). 

[36] Z. X. Shen et al. Phys. Rev. Letters 70, 1553 (1993), and papers cited therein; ibid. 
private communication. 

[37] M. J. Rice, and Y. R. Wang, preprint. 

[38] N. Nagaosa, and P. A. Lee, Phys. Rev. Lett. 64, 2450 (1990). 
[39] V. Kalmeyer, and B. Laughlin. Phys. Rev. Lett. 59, 2095 (1987). 
[40] P .B. Littlewood, Phys. Rev. B 42, 10075 (1990). 

[41] R. Micnas, J. Ranninger, and S. Robaszkiewicz, Rev. Mod. Phys. 62, 113 (1990). 
[42] D. van der Marel, and J. E. Mooij, Phys. Rev. B. 45, 9940 (1992). 

16 



FIGURES 

FIG. 1. First Brillouin-zone of a square lattice, with the occupied states indicated as the shaded 
area. The lozenge indicates the perfectly nested Fermi surface. 

FIG. 2. Density of states at the Fermi energy in units of 1/W as a function of electron occu- 
pation number 

FIG. 3. Free energy difference with the normal state of the s*-wave (solid) and <i-wave(long 
dashed) superconducting state and of the anti-ferromagnetic state (short dashed curve) with 
Jo = 0.6W and if = 0. Energies are in units of W. 

FIG. 4. Phase diagram in the Jq-u plane, where n is the number of electrons per unit cell. 

FIG. 5. Free energy difference with the normal state of the s*-wave (solid) and <i-wave(long 
dashed) superconducting state and of the anti-ferromagnetic state (short dashed curve) with 
Jo = 0.6 W and if = — 0.7t. Energies are in units of W. 

FIG. 6. Solid curve: T c /Jq calculated for the d-wave channel of the exchange-only model with 
if = and 1 electron per site. The same curve is obtained for if = 0.7 with 0.7 electron per site. 
Open lozenges: T c of the s*-wave channel with the latter parameters. Dotted curve: T C /\U\ versus 
\U\/W for the negative U Hubbard model taking a square DOS. 

FIG. 7. Phase diagram in the temperature-density plane with the parameters Jo = O.leV, 
W = 8t = O.beV, and t'/t = —0.7. The curves are interrupted in the part near the middle, where 
the present analysis is physically meaningless. The mean field Neel temperature at half filling, 
using the same value for Jq as in the metallic regime, is indicated as a clover. 
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